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RAs, Beet, RE! r3 

《1. 中 国 科学 院 国家 天 文 台 ， 北 京 100101; 2. 中 国 科学 院 大 学 ， 北 京 100049; 3. 中 国 

科学 院 月 球 与 深 空 探 测 重点 实验 室 , 北京 100101; 4. 清 华 大 学 工程 物理 系 , 北京 100084) 
W 要 : HAO TIAN E MAEAEA RER, AARAM EARI. Rin, HR 
四 号 平台 存在 约 10-5W/(m2*Hz) 级 别 的 强 干 扰 ， 并 且 干 扰 在 每 道 时 域 数据 中 都 存在 明显 的 差 
Jt, 这 大 大 前 弱 了 低频 射电 频谱 仪 的 天 文 观测 灵敏 度 。 为 此 ， 本 论文 从 两 组 信号 的 相关 性 出 
发 ， 提 出 基于 CLEAN 算法 ， 借 助 互 相关 功率 谱 、 传 立 叶 级 数 展开 等 工具 ， 把 低频 射电 频谱 
BLA, B. C 天 线 的 时 域 观测 数据 切 分 为 强 相关 的 CLEAN 模型 信号 和 部 分 相关 的 残余 信号 
两 个 部 分 。 其 中 CLEAN 模型 信号 主要 由 平台 干扰 信号 和 可 能 的 低频 强 射电 爆发 组 成 ; 残余 
信号 由 接收 机 噪声 、 未 扣除 干净 的 平台 干扰 信号 和 常规 的 低频 射电 信号 组 成 。 将 该 方法 应 用 
到 实际 数据 中 ， 结 果 表 明 嫦 娥 四 号 低频 射电 频谱 仪 的 未 积分 灵敏 度 可 以 提高 约 8 个 数量 级 ， 
达到 103W/(m2*Hz) 水 平 。 在 此 基础 上 ， 基于 对 平台 干扰 信号 中 确定 性 成 分 和 宽带 随机 成 分 
的 分 别处 理 ， 再 借助 低频 射电 爆发 信号 和 平台 干扰 信号 在 功率 谱 上 不 同 的 表现 ， 以 及 常规 低 
频 射 电 天 文 信 号 受 月 球 自转 调制 等 信息 ,将 来 科学 分 析 工 作 的 重点 是 进一步 处 理 CLEAN 模 
型 信号 和 残余 信号 ， 以 期 发 现 低频 强 射 电 天 文 爆发 信号 ， 乃 至 对 全 天 区 进行 粗略 的 成 像 。 
关键 词 : 嫦娥 四 号 低频 射电 频谱 仪 ， 平 台 干 扰 分 离 ，CLEAN 算法 。 
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1 引言 

嫦娥 四 号 低频 射电 频谱 仪 (Low frequency radio spectrum analyzer 简称 LFRS) 放置 
在 月 球 背 面 。 由 于 月 球 屏 蔽 了 大 量 来 自 地 球 的 无 线 电 干 扰 ， 因 此 LFRS 具有 得 天 独 厚 的 观测 
条 件 ， 可 以 对 太阳 等 宇宙 天 体 的 射电 信号 进行 观测 。" 但是， 由 于 着 陆 器 平台 干扰 信号 比 接 
收 到 的 低频 射电 天 文 信号 大 很 多 ,如果 不 分 离 平 台 干 扰 信 号 , 将 无 法 得 到 有 效 的 低频 射电 天 
文 信号 数据 。 因 此 ， 需 要 提出 一 种 有 效 分 离 平 台 干 扰 信 号 的 方法 ， 从 而 实现 对 干扰 的 抑制 。 
在 射电 天 文 领域 ， 传 统 抑 制 干扰 的 方法 有 很 多 ,包括 对 空白 脉冲 干扰 的 时 间 限 制 ”"， 自 
适应 频率 选择 滤波 "， 参 数 信号 相 减 ”"， 自 适应 波束 成 形 “”， 时 间 自 适应 滤波 "， 带 有 参考 天 
线 的 频 域 后 相关 处 理 ”"， 子 空间 投影 的 空间 滤波 "等 等 。 到 目前 为 止 ， 运用 比较 广泛 的 算法 
是 单 通 道 总 功率 变化 检测 器 ， 它 使 得 相关 器 的 输出 得 到 减 小 。Friedman (1996) 基于 对 功率 
变化 的 检测 ， 在 RATAN600 上 实现 了 改进 的 功率 检测 器 ”。Weber 5& A. (1997 F) 提出 在 所 
滞后 的 情形 下 使 用 量化 相关 性 来 测试 干扰 的 存在 ""。Maslakovic 等 人 1996 年 提出 了 一 种 
基于 小 波 分 解 的 检测 器 ”…。 这 些 都 是 单 通道 检测 器 , 它们 不 利用 干扰 的 空间 特性 分 析 。.Kasper， 
Chute 和 Routledge 在 1982 年 提出 了 将 多 个 望远镜 结合 使 用 以 改善 探测 质量 的 方法 ， 此 方 
法 用 于 低频 干涉 测量 , 在 这 之 中 还 提出 了 一 种 基于 交叉 频谱 时 间 行 为 的 鲁 棒 数 据 检 查 方法 ”。 
Barnbaum&Bradley(1998 ) 考 虑 了 自 适 应 过 滤 技 术 , 他们 提议 使 用 参考 天 线 和 最 小 均 方 (LMS) 
类 型 算法 消除 Green-Bank 射电 望远镜 的 干扰 "”。W. A. BAAN 提出 了 一 种 针对 韦 斯 特 博克 综合 
射电 望远镜 的 射电 干扰 消除 子 系统 "。 这 些 都 是 在 射电 干扰 消除 方面 典型 的 例子 。 

上 面 提 到 的 射电 干扰 信号 消除 算法 主要 分 为 三 类 ; 第 一 类 试图 去 通过 线性 的 方式 来 区 分 
射电 干扰 信号 的 特点 ， 比 如 奇异 值 分 解 SVD) UU, SATA PT (PCA) "等 等 。 如 果 射 电 干 
扰 在 时 间 和 频率 上 表现 出 重复 的 模式 ,这些 方法 效果 很 好 ,但 无 法 处 理 更 多 随机 信号 , 例如 
由 卫星 引起 的 不 规则 信号 。 第 二 类 是 基于 闵 值 的 算法 , 例如 CUMSUM™ 和 SUMTHRESHOLD", J 
管 它 们 很 简单 , 但 是 这 些 方法 相当 可 靠 并 且 非 常 有 效 。 第 三 类 使 用 传统 的 监督 机 器 学 习 技 术 
(例如 kK 最 近邻 和 高 斯 混合 模型 ) 对 射电 干扰 信号 进行 聚 类 。 

嫦娥 四 号 LFRS 所 收 到 的 平台 干扰 信号 很 强 ,和 太阳 爆发 峰值 的 量 级 10”W/ (Hz) 相当， 


202105.00045v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


比 普通 的 太阳 爆发 信号 10 W/ mH 大 三 个 量 级 ， 比 太阳 宁静 时 的 信号 10 7W/ (mHz) 
(1S.F.U) 大 七 个 量 级 。 之 前 仪器 设计 时 推荐 使 用 的 谱 减 法 ""， 虽 然 可 以 将 信号 的 信 噪 比 提 
升 约 30dB， 但 是 仍然 不 足以 探测 有 效 的 射电 天 文 信 号 。 为 此 ， 本 文 提出 了 基于 CLEAN 算法 
的 平台 干扰 信号 分 离 方 法 ， 利 用 LFR 三 根 天 线 同时 观测 信号 存在 相关 性 的 特点 ， 并 结合 
线 之 间 数 据 的 互相 关 功 率 谱 、 傅 立 叶 级 数 展开 等 工具 ， 有 效 地 分 离 平 台 干 扰 信号 ， 从 而 实现 
对 干扰 的 抑制 。 
本 文 的 第 2 小 节 主 要 介绍 LFRS 的 工作 原理 及 接收 信号 构成 成 分 ， 第 3 小 节 介 绍 LFRS 
平台 干扰 信号 分 离 的 方法 及 算法 流程 ， 第 4 小节 展示 新 算法 对 LFRS 实际 数据 的 处 理 结果 及 
效果 分 析 ， 第 5 小 节 是 讨论 ， 最 后 的 第 6 小 节 是 结论 。 
2 LFRS 的 工作 原理 及 接收 信号 的 构成 
2. 1LFRS 简要 介绍 
太阳 和 宇宙 中 的 其 它 来 波 信号 是 一 个 矢量 信号 。 根据 矢 量 场 的 理论 , 来 波 信号 可 以 分 解 
为 任意 三 个 互相 垂直 方向 的 分 量 ,， 只 要 将 这 三 个 分 量 接收 下 来 ,进行 合成 ,就 可 以 得 到 矢量 
的 大 小 和 方向 。 通 过 后 期 处 理 ， 可 以 获得 电场 的 频谱 和 时 变 信息 。 
LFRS 的 天 线 分 布 如 图 1 所 示 : 三 分 量 有 源 天 线 接收 宇宙 的 来 波 电场 。 三 根 天 线 分 别 标 
注 为 “天 线 A”、“ 天 线 B” 和 “天 线 C”， 两 两 正 交 ， 还 有 一 根 短 的 “天 线 D” 主要 用 于 接收 
平台 干扰 信号 。 为 天 线 数据 处 理 做 一 个 参考 。 


图 1 着 陆 器 低频 射电 频谱 仪 的 天 线 位置 分 布 
Fig. 1 Antenna location distribution of lander low-frequency radio spectrum analyzer 


LFRS 具有 高 频 和 低频 两 个 频率 接收 波段 , 高 频 对 应 的 频率 范围 是 1-40MHz, 低频 对 应 的 


频率 范围 是 0. 1-2MHz。 接 收 机 灵敏 度 为 6nV/VHz， 高 频 的 频率 分 辨 率 为 100kHz， 低 频 的 频 
率 分 辩 率 为 5kHz。 


LFRS 主要 由 四 根 接收 天 线 、 前 置 放大 器 、 电 子 学 单元 、 电 缆 组 件 等 部 分 组 成 。 有 源 天 
线 安装 于 卫星 本 体外 ， 其 余 组 件 安装 于 卫星 本 体内 。 电 子 学 单元 由 控制 器 、 配 电器 、 基 准时 
钟 模块 、 四 通道 接收 机 、 内 定 标 模块 、 存 储 单元 、 通 讯 接口 、 数 传 接口 等 组 成 。LFRS 内 部 
各 单元 之 间 及 与 数 管 分 系统 、 供 配 电 分 系统 之 间 的 接口 关系 如 图 2 所 示 。 
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关系 示意 图 


Fig. 2 Schematic diagram of low-frequency radio spectrum analyzer, wit 
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LFRS 有 六 种 工作 模式 ， 分 别 为 1) 待机 工作 模式 2) 内 定 标 工作 模式 3) 时 频 对 比 工作 
模式 4) 频谱 巡查 工作 模式 5) 三 天 线 工 作 模 式 6) 时 域 工作 模式 。 其 中 运用 最 广泛 的 是 三 


天 线 工 作 模式 ， 在 此 模式 下 ，LFRS 分 别 采 集 四 根 天 线 
处 理 的 数据 均 对 应 三 天 线 工作 模式 。 


四 个 通道 的 科学 数据 ， 本 文 进行 读 取 


LFRS 的 科学 数据 分 为 0A 级 、0B 级 、01 级 、2A 级 、2B 级 、2C 级 共 六 级 。 其 中 预 处 理 数 


据 经 过 分 路 解 帧 和 分 包 ， 生 成 0A 级 数据 ;0A 级 数据 经 过 排序 、 去 习 


E 复 、 两 站 优化 拼接 、 去 


源 包 包头 、 根 据 采样 周期 形成 LFRS 数据 块 ， 生 成 OB 级 数据 产品 ， 多 个 OB 级 数据 文件 拼合 
成 一 个 探测 周期 ， 生 成 01 级 数据 产品 ， 对 01 级 数据 进行 傅立叶 变化 、 定 标 等 操作 ， 生 成 
2A 级 数据 产品 ， 对 2A 级 数据 加 入 太阳 位 置 等 信息 ， 生 成 2B 级 数据 产品 ， 对 2B 级 数据 进行 
对 消 噪声 处 理 ， 生 成 2C 级 数据 产品 。 其 中 1 级 数据 为 时 域 数 据 ，2 级 数据 为 频 域 数据 。 本 


论文 处 理 所 用 的 原始 数据 均 为 1 级 数据 。 
2. 2 接收 信号 的 构成 成 分 


LFRS 收 到 的 信号 主要 由 三 部 分 构成 ， 一 部 分 是 平台 干扰 信号 ， 


最 后 是 低频 射电 天 文 信号 。 


部 分 是 接收 机 噪声 ， 


平台 干扰 信号 的 特点 是 幅度 相对 较 强 , 两 两 天 线 之 间 的 相关 性 高 。 如 果 干 扰 源 位 置 确定 ， 


各 天 线 接收 信号 之 间 的 强度 比例 不 随时 间 变 化 。 且 来 自 电 子 设备 的 干扰 信号 , 一 般 可 分 为 固 
部 分 是 宽带 的 随机 噪声 信号 ， 这 是 由 
于 EMC 兼容 问题 带 来 的 。 最 终 的 平台 于 扰 信 号 如 图 6 所 示 ， 既 有 确 


pet 


定 频率 的 发 射 信号 ， 如 通信 发 射 机 ， 时 钟 信号 等 ， 另 


带 随 机 性 噪声 成 分 。 


oH 


E 变 化 的 成 分 ， 也 有 宽 


低频 射电 天 文 信号 也 具有 一 定 的 相关 性 ， 除 此 之 外 ,还 有 一 些 天 体 源 所 特有 的 性 质 。 比 
W, 低频 太阳 射电 信号 其 共有 了 瞬时 性 和 频谱 漂移 等 特点 , 为 此 其 具有 很 强 的 频谱 可 识别 特征 
"9. Ql 3 所 示 。 再 比如 ， 快 速射 电 爆 (Fast Radio Burst, FRB) FRB200428 的 持续 时 间 
其 辐射 强度 足以 被 LFRS 探测 到 。 另 
外 ， 由 于 仪器 指向 随 月 球 的 自转 而 改变 ， 低 频 射 电 天 文 信号 在 LFRS 各 天 线 之 间 的 强度 比例 


约 0. 6lms， 能 量 通 量 (Fluence) 达到 1.5 MJy ms”, 


随时 间 呈 现 有 规律 的 变化 。 所 有 这 些 特性 ， 都 有 助 在 绰 
射电 天 文 信号 。 
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3: 低 频 太 阳 射 电信 号 频谱 


Fig.3 Low frequency spectrum of solar radio signal 
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接收 机 噪声 的 特点 是 天 线 两 两 之 间 无 相关 性 。 可 以 利用 该 特性 ,通过 计算 天 线 之 间 数 据 


的 相关 性 来 评估 平台 干扰 信号 的 分 离 程度 。 
图 4 展示 了 LFRSB 天 线 观测 到 的 一 个 频谱 实例 。 由 图 可 以 得 知 ， 


频谱 在 1MHz 频率 附近 


的 流量 密度 大 约 为 10-35W/ (mx*Hz)， 与 太阳 爆发 峰值 流量 密度 10-15W/ (mHz) 相当 ， 比 太阳 
爆发 流量 密度 10- 8W/ (mHz) X f 2-3 个 量 级 , 说 明 着 陆 器 平台 干扰 比 我 们 希望 探测 的 太阳 
爆发 信号 要 大 许多 。 如 果 不 通过 一 定 的 方法 分 离 平 台 干 扰 信 号 , 将 无 法 有 效 地 提取 出 想 要 的 


低频 射电 天 文 信号 。 


The Flux Density of antena B before CLEAN 
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图 4 B 天 线 采 集 到 的 高 频 波 段 的 一 个 频谱 实例 ， 对 应 第 23 A 


二 道 数据 


Fig. 4 An example of the spectrum of antennaB at the high-frequency band, corres 


data on the 23rd moon day. 


3 LFRS 平台 干扰 信和 号 的 分 离 方法 


ponding to the second channel 


低频 射电 源 与 平台 干扰 源 会 在 LFRS 的 3 个 天 线 上 激发 出 互相 关联 的 接收 信号 。 本 论文 
的 目的 是 从 两 组 信号 的 相关 性 出 发 ， 借 助 CLEAN 方法 ， 以 及 互相 关 功 率 谱 、 傅 立 叶 级 数 展开 


等 工具 ， 把 LFRS 观测 信号 中 的 强 相关 成 分 与 部 分 相关 成 分 进行 分 离 ， 


以 便 进行 后 续 的 科学 
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分 析 。 分 离 之 后 ， 原 始 的 LERS 时 域 信号 被 切 分 为 两 个 成 分 ， 即 强 相关 的 CLEAN 模型 信号 
部 分 相关 的 残余 信号 。 在 这 个 过 程 中 ， 原 始 信号 的 信息 没有 丢失 或 损伤 。 
CLEAN 模型 信号 主要 由 平台 干扰 信号 组 成 。 当 低频 强 射电 爆发 信号 到 来 时 ， 也 会 进入 
CLEAN 模型 信号 。 由 于 数 传 的 限制 ，LFRS 大 约 以 1 秒 钟 为 间隔 采集 一 道 时 域 数 据 ,， 采 集 间 隔 
并 不 严格 相等 。 每 道 数据 由 4096 个 数据 组 成 ， 采 样 频率 为 100MHz， 采 样 持续 时 间 只 有 大 约 
40 微 秒 。 因 此 ， 不 同道 时 域 数据 有 明显 的 差别 (参见 图 6 与 图 7)。 
为 了 寻找 有 效 的 低频 强 射电 爆发 信号 , 需要 提取 每 一 道 时 域 数据 的 CLEAN 模型 信号 , 并 
对 其 中 的 确定 性 成 分 和 宽 带 随 机 成 分 分 别 进行 整理 。 例如 ,对 于 CLEAN 模型 信号 中 的 确定 性 
时 域 变 化 结构 或 其 频谱 成 分 ， 可 以 把 它们 当 作 某 种 系统 误差 ， 同 样 借助 CLEAN 算法 , 将 其 扣 
除 。 剩 下 的 CLEAN 模型 信号 ， 由 宽带 随机 信号 组 成 ， 可 以 通过 平均 多 道 数 据 的 方式 ， 降 低 其 
涨 落水 平 ， 提 高 探测 低频 强 射 电 天 文 爆发 的 灵敏 度 
残余 信号 主要 由 不 相关 的 接收 机 噪声 、 未 分 离 出 去 的 平台 干扰 信号 、 以 及 常规 低频 射电 
天 文 信号 组 成 。 残 余 信号 的 特征 类 似 于 随机 白 噪声 ， 里 面 没 有 确定 性 的 信号 成 分 ， 可 以 借助 
于 数理 统计 工具 来 进行 分 析 。 例如， 可 以 借鉴 微波 辐射 计 的 思想 , 通过 积分 平均 的 方式 进 一 
` 步 提 高 残余 信号 的 探测 灵敏 度 。 
~ 我 们 的 最 终 目标 是 用 嫦娥 四 号 LFRS 的 观测 数据 来 搜索 低频 射电 天 文 信号 ， 并 进行 深入 
t 的 科学 研究 。 此 时 ， 无 论 是 CLEAN 模型 信号 还 是 残余 信号 都 是 缺 一 不 可 的 。 
时 标 短 (分 钟 到 天 量 级 )、 幅 度 强 (MJy 量 级 ) 的 低频 强 射电 爆发 往往 隐藏 在 CLEAN 模 
型 信号 中 。 对 这 类 射电 爆发 信号 的 提取 离 不 开 宁 静 状 态 的 CLEAN 模型 信号 的 分 类 整理 。 特别 
强 的 低频 射电 天 文 爆发 ， 在 单 道 数据 中 就 有 明显 异常 的 表现 。 此 时 ， 可 以 借助 交叉 相关 ， 直 
接 扣 除 该 道 数据 对 应 的 宁静 状态 的 CLEAN 模型 信号 ， 获 得 射电 爆发 信号 的 时 域 、 频 谱 特 征 。 
如 果 射 电 爆 发 强度 与 平台 干扰 中 的 宽带 随机 成 分 接近 , 那么 先 扣 除 CLEAN 模型 信号 中 的 确定 
性 成 分 ， 再 通过 多 道 数据 积分 平均 的 方式 ,提高 探测 灵敏 度 ， 最终 获 得 这 类 低频 射电 天 文 爆 
发 信号 的 时 域 、 频 谱 特 征 。 
一 些 时 标 长 (超过 月 量 级 )、 幅 度 弱 的 低频 射电 天 文 信号 则 隐藏 在 残余 信号 中 。 为 了 探 
测 这 类 弱 信 号 ， 首 先 要 借助 星 载 定 标 装置 或 其 它 外 部 稳定 的 定 标 源 ， 确 认 接 收 机 工作 稳定 。 
之 后 ,可 以 直接 从 残余 信号 出 发 ， 同样 借 助 微波 辐射 计 的 思路 ， 通 过 积分 平均 ， 提 高 残余 信 
号 的 探测 灵敏 度 。 理 论 上 讲 ，LFRS 三 个 天 线 残 余 信号 的 积分 值 会 随 着 月 球 自转 而 发 生 调制 
变化 。 这 样 的 调制 信息 甚至 可 以 用 来 对 全 天 进行 粗略 的 成 像 。 
3. 1 用 CLEAN 算法 分 离 平 台 干 扰 信 号 的 基本 思路 
设 LFRS 中 一 个 天 线 收 到 的 信号 为 S(D，, 它 包含 了 平台 干扰 信号 Kb, 字 宙 身 
以 及 接收 机 噪声 N(t)。 这 些 信号 之 间 的 关系 如 下 : 
Sa(t) = aaO) + PADCE) + Nal) 
Sg(t) = ag (HICE) + Bg(t)C(t) + Ne(t) (1) 
Sc(t) = ac (ÐI) + POCE) +N) 
其 中 Sj(t)、Sp(t)、Sc(t) 分 别 为 天 线 A、B、C 接收 到 的 信号 ; NO Na(t). Nc(t)2) 
为 天 线 A、B、C 对 应 接收 机 的 本 地 噪声 ， 它 们 之 间 相 互 独立 ;Qay(t)、ap(t)、ac(t) 为 平台 
扰 信号 在 天 线 A、B、C 上 的 投射 系数 ，P4(t)、 "app toc pea C 
上 的 投射 系数 。 
在 2.2 小 节 中 提 到 ,平台 干扰 信号 比 低频 射电 天 文 信号 大 约 3 一 10 eR. Alt, 2B 
台 干 扰 是 进行 有 效 天 文 观测 的 一 个 基本 前 提 。 
经 典 的 CLEAN 算法 在 脏 图 或 残余 图 搜索 最 大 值 ,此 位 置 处 包含 部 分 真实 目标 , 用 减 点 源 
扩散 函数 (PSF) 的 方式 将 这 部 分 真实 目标 提取 出 来 ， 如 此 不 断 和 迭代 ， 直 到 残余 图 中 只 剩 噪 
结构 为 止 ”。CLEAN 算法 同样 可 以 应 用 到 时 域 信号 上 ， 比 如 利用 CLEAN 方法 对 时 间 序 列 进 
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个 加 的 互相 关 功 率 谱 包 含 了 平台 干扰 信号 对 应 的 频谱 成 分 , 同时 夺 括 了 振幅 和 相位 信息 。 
根据 这 些 信息 , 我 们 借助 CLEAN 算法 和 傅立叶 级 数 展开 公式 , 可 以 逐一 将 某 个 频率 成 分 对 应 
的 三 角 函 数 信号 从 时 域 数据 中 扣除 。 考 虑 到 每 个 天 线 及 其 接收 机 的 增益 并 不 一 致 , 在 扣除 时 


每 个 天 线 的 CLEAN 因子 是 不 一 样 的 。 
坦 立 叶 变化 所 对 应 的 时 域 及 频 域 
会 在 傅立叶 变换 时 出 现 频 谱 泄 漏 ， 导 致 频 率 估计 出 现 偏差 。 
满足 首尾 连续 的 条 件 ， 正 如 图 6 
函数 的 方式 ， 来 避免 频谱 泄漏 。 但 是 加 窗 函 数 后 会 导致 信号 发 生 扭 曲 ， 幅 度 值 发 生变 化 ， 
立 叶 反 变 换 的 方式 来 进行 信号 分 离 ， 则 会 导致 较 大 的 偏差 。 为 此 ,我 们 
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言 号 的 信 噪 比 ， 并 利用 信号 随时 间 的 变化 特性 来 进 
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天 文 信号 ， 从 而 实现 对 干扰 的 抑制 。 


3. 2 CLEAN 算法 的 流程 


按照 3. 1 子 小 节 的 思路 , 我 们 给 出 用 CLEAN & 
明 清晰 的 展示 算法 ,我 们 先 考 虑 两 根 天 线 (分别 记 
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昌 度 往往 小 于 噪声 水 平 ， 加 上 小 部 分 
， 可 以 通过 对 时 域 信号 积分 平均 的 手段 ， 提 高 残余 
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其 中 M4 和 ME5 分 别 为 FM(w) 与 PE(w) 的 平均 振幅 ，M4 人 为 P48(w) 的 平均 振幅 。 
FEO = fA) — EMi, cos (Wnt + Pm) 


FECE) = fE(t) — 6M2 cos (wnt + Pm) (3) 
fioa = fioa) 十 6M4 cos (ct + Pm) 
Foa) = fo (D + 6 Mp cos (Wnt + Pm) (4) 


AK SUE EF] SB IEDR MI MM, Be EPR EF (0) RIF P (eo) FEST cos, KEATS 
其 原因 是 当 信 品 比较 低 时 , 频谱 的 振幅 由 噪声 谱 振 幅 主 导 , Tft EL ROUGH PPP (oo) FER Om, 
处 的 振幅 M 铝 受 噪声 的 影响 较 小 。 
在 这 里 ， 我 们 假定 天 线 A 和 B 的 时 域 信 号 天 (9D 和 产 (9 是 同步 的 。 若 由 于 线路 补偿 没有 
做 好 ， 两 个 时 域 信 号 之 间 存 在 时 间 延 迟 ， 那 么 需要 通过 交叉 相关 等 方法 计算 该 时 间 延 迟 ， 并 
对 时 域 信号 做 延迟 修正 ， 然 后 再 用 CLEAN 算法 做 平台 干扰 信号 的 提取 。 

通过 CLEAN 算法 分 离 平台 干扰 信号 之 后 ,得 到 的 残余 信号 主要 成 分 是 接收 机 噪声 ， 以 及 
隐藏 其 中 的 宇宙 射电 信号 , 还 有 少 部 分 没有 扣除 干净 的 平台 干扰 信号 。 此 时 两 个 天 线 残 余 信 
号 之 间 的 相关 系数 比较 低 ， 例 如 绝对 值 小 于 0.3。 这 样 ， 通 过 计算 信 噪 比 发 现 ， 残 余 信 号 
的 平台 干扰 信号 和 接收 机 噪声 基本 在 同一 个 量 级 水 平 上 。 再 往 后 , 可 以 通过 积分 平均 的 方式 
来 提高 残余 信号 的 信 噪 比 , 这 与 传统 辐射 计 的 工作 思路 一 致 。 至 于 平台 干扰 信号 和 宇宙 射电 
信号 的 区 分 ， 则 要 借助 三 个 天 线 信 和 号 随 频谱 、 时 间 和 来 波 方向 变化 的 规律 。 


4 LFRS 实例 数据 处 理 
4. 1LFRS 数据 处 理 流程 

LFRS 原始 数据 由 嫦娥 四 号 地 面 应 用 系统 数据 发 布 中 心 提供 。 数 据 放 在 嫦娥 四 号 低频 射 
电 频 谱 仪 载荷 相关 文件 夹 下 ， 按 照 0A、0B、01、2A、2B 进行 分 级 排列 ， 每 一 级 数据 均 按 照 
每 一 个 月 昼 一 个 文件 夹 进行 排列 ， 每 一 个 月 恒通 常 有 14 个 文件 ， 其 中 7 个 是 数据 文件 ，7 
个 是 数据 采集 参数 文件 。7 个 数据 文件 对 应 的 是 在 一 个 月 昼 中 7 次 开机 时 间 段 所 采集 到 的 科 
学 数据 ， 我 们 在 研究 不 同时 间 段 的 数据 时 ， 分 别 对 这 7 个 文件 进行 读 取 ， 并 结合 相应 的 7 
个 数据 采集 参数 文件 ， 对 每 一 个 时 段 的 采集 参数 进行 确定 。 

以 第 23 月 尽 高 频 1 级 时 域 数 据 的 第 一 道 数 据 为 例 。 对 应 的 数据 文件 名 为 
CE4 GRAS LFRS-TR SCI P 20201014020000 20201014140000 0135 B. 01， 其 中 GRAS 代表 地 
看 应 用 系统 ，SCI 代表 科学 数据 ，0135 代表 这 个 文件 是 总 共 第 135 个 数据 文件 ， 数 据 结尾 的 
01 扩展 名 代表 数据 是 1 级 数据 , 为 时 域 数据 。 从 文件 名 可 以 看 出 采集 的 数据 时 间 段 是 从 2020 
年 10 月 14 日 2 点 0 分 0 秒 到 2020 年 10 月 14 日 14 点 0 分 0 秒 , 这 个 时 间 是 世界 时 , 转换 
为 北京 时 间 还 需要 加 8 小 时 进行 换算 。 
在 数据 文件 中 ， 数 据 按照 道 的 顺序 依次 排列 ， 每 一 道 数据 由 数据 头 〈79 字 节 ) 和 实际 
科学 数据 (8192 字 节 ) 组 成 。 其 中 数据 头 包 含 了 仪器 的 相关 参数 ， 科 学 数据 共有 4096 个 科 
学 数据 点 ， 每 个 数据 2 个 字 节 ， 由 16 位 无 符号 整数 组 成 ， 高 字 节 在 前 ， 低 字 节 在 后 。 读 取 
后 还 需要 将 整数 电压 转化 为 浮 点 数 电 压 值 , 具体 数据 形式 可 参考 文档 《 探 月 工程 剥 娥 四 号 任 
务 着 陆 器 低频 射电 频谱 仪 数据 预 处 理 方法 设计 》"“。 第 一 道 数据 一 共 采 集 了 4096 个 点 ， 高 
频数 据 的 采样 间隔 为 0.01hs， 低 频数 据 的 采样 间隔 为 0.2hs。 

我 们 是 通过 python 编程 的 方式 ， 将 数据 按照 道 的 顺序 依次 读 取出 来 ， 将 每 一 道 4096 
个 点 排 成 一 行 ， 生 成 一 个 科学 数据 和 矩阵。 处 理 的 时 候 ， 提 取 每 一 行 〈 每 一 道 ) 数据 ， 分 别 用 
第 3. 2 子 小 节 所 述 的 CLEAN 方法 进行 处 理 。 
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4. 2LFRS 部 分 处 理 结果 展示 

在 本 子 小 节 ， 我 们 将 分 别 展 示 CLEAN 前 后 天 线 时 域 数据 (图 6、 图 11)， 频 域 数据 (图 
8, El 130, 流量 密度 随 频 率 分 布 的 曲线 (图 9)，CLEAN 后 model 信号 (图 12)， 互 相关 功率 
WE (图 10、 图 14)， 相 关系 数 随 迭 代 次 数 变化 曲线 〈 图 15)， 并 对 相应 的 结果 进行 说 明 。 
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图 6 A、B、C、D 天 线 CLEAN 前 的 第 23 月 昼 第 一 道 原始 时 域 数 据 
Fig.6 The first time-domain data, before CLEAN, of antenna A, B, C, and D on the 23rd moon day. (a)Time domain data 


of antenna A; (b) Time domain data of antenna B; (c) Time domain data of antenna C; (d) Time domain data of antenna 
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The Original Data of C Antenna The Original Data of D Antenna 
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17 A、B、C、D 天 线 CLEAN 前 的 第 23 月 屋 第 二 道 原始 时 域 数据 


Fig.7 The second time domain data, before CLEAN, of antenna A，B，C，and D on the 23rd moon day. (a)Time domain 
data of antenna A; (b) Time domain data of antenna B; (c) Time domain data of antenna C; (d) Time domain data of 
antenna D 

从 CLEAN 前 的 时 域 数 据 (图 6) 中 可 以 看 到 ，A 天 线 和 B 天 线 的 信号 值 范围 基本 接近 ， 
均 为 -0. 5V~0. 5V 左右 。C 天 线 数据 比 A、B 天 线 略 大 。 三 个 天 线 的 信号 随时 间 变 化 的 结构 具 
了 较 高 的 相似 性 ，A、B 之 间 的 相关 系数 高 达 0.965，B、C 之 间 的 相关 系数 为 0. 823, A C 
之 间 的 相关 系数 为 0.915。 同 时 ， 数 据 中 存在 的 来 回 震 荡 结 构 主 要 成 分 不 是 接收 机 噪声 ， 
为 接收 机 噪声 之 间 无 相关 性 ,A、B、C 天 线 信 号 的 均值 分 别 为 -0. 0149V、-0. 0143V 和 -0. 0032V, 
均 接 近 于 0。 它 们 的 标准 差 c 分 别 为 0. 189V、0. 138V 和 0. 409V， 存 在 明显 的 差别 。 这 种 差 
别 可 能 是 由 于 三 天 线 的 增益 不 同 ,或 者 平台 干扰 在 三 天 线 上 的 投射 系数 不 一 致 , 具体 原因 还 
了 赖 于 后 续 细 致 的 数据 分 析 。 从 CLEAN 前 第 二 道 时 域 数据 (图 7) 中 可 以 看 到 ， 第 一 道 数据 
与 第 二 道 数据 时 域 信号 有 较为 明显 的 差别 。 

从 四 根 天 线 的 数据 来 看 ， 发 现 D 天 线 的 响应 值 比 其 余 三 根 天 线 小 很 多 ,这 是 由 于 DD 天线 
长 度 只 有 20cm， 比 A、B、C 长 度 (5m) 小 很 多 造成 的 。 同 时 ，D 天 线 信号 与 其 余 三 根 天 线 信 
号 在 形状 上 也 有 具有 明显 的 差别 , 说明 D 天 线 的 频谱 响应 明显 异 于 其 余 三 根 天 线 , 并 且 这 种 差 
异 是 非 线 性 的 。 从 天 线 之 间 的 相关 系数 也 可 以 看 出 ，A、B、C 三 根 天 线 与 D 天 线 的 数据 相关 
系数 分 别 为 0. 440，0. 397 和 0. 298， 显 著 低 于 A、B、C 天 线 信 号 之 间 的 相关 系数 。 因 此 ， 
我 们 放弃 使 用 D 天 线 数据 来 进行 平台 干扰 信号 的 分 离 。 


The Spectrum of B Antenna Before CLEAN The Spectrum of A Antenna Before CLEAN 
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Ca) A 天 线 时 域 数据 CLEAN 前 的 频谱 Cb) B 天 线 时 域 数据 CLEAN 前 的 频谱 


图 8 A、B 天 线 CLEAN 前 第 23 月 昼 第 一 道 数据 的 频谱 


Fig.8 The spectrum of the first data, before CLEAN, of antenna A and B on the 23rd moon day. (a) The spectrum, before 


CLEAN, of antenna A; (b) The spectrum, before CLEAN, of antenna B; 
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The Flux Density of antena A before CLEAN The Flux Density of antena B before CLEAN 
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Ca) A 天 线 CLEAN 前 流量 密度 (b) B 天 线 CLEAN 前 流量 密度 


图 9 A、B 天 线 第 23 月 昼 第 一 道 数据 CLEAN 前 的 流量 密度 


Fig.9 The first flux density data, before CLEAN, of antennas A and B on the 23rd moon day. (a)Flux density of antenna 


A before CLEAN; (b) Flux density of antenna B before CLEAN 
从 两 天 线 CLEAN 前 的 频谱 (图 8) 可 以 看 到 ， 频 谱 的 幅度 谱 强度 分 布 在 10-1VY~10-Y 之 
间 ， 低 频 成 分 强度 较 大 ， 并 且 频 谱 中 呈现 出 复杂 的 结构 。 将 此 时 的 电压 频谱 转化 为 流量 密度 
频谱 (图 9) 后 ,发现 两 天 线 较 低频 段 的 流量 密度 峰值 可 超过 10-16W/ (m2 Hz)， 比 太阳 爆 
发 信号 的 流量 密度 10-18W/(m?: Hz) 还 要 大 2 个 量 级 。 说 明 直 接 读 取 的 时 域 信号 如 果 不 经 过 
去 除 平台 干扰 的 处 理 ， 我 们 很 难 发 现 宇宙 射电 爆发 信号 。 


The Cross-correlation spectrum Before CLEAN 


T 
1.0 10.75 20.5 30.25 40.0 
frequency(MHz) 


图 10 A. B 天 线 CLEAN 前 的 互相 关 功 率 谱 
Fig. 10 Cross power Spectrum of antennas A and B before CLEAN 
两 天 线 在 CLEAN 前 的 互相 关 功 率 谱 〈 图 10) 显示 其 中 存在 一 些 较为 复杂 的 结构 。 此 外 ， 
与 单 天 线 频谱 相 比较 ， 两 天 线 互相 关 功 率 谱 的 信 噪 比 更 高 。 对 比 图 10 和 图 8， 可 以 发 现 单 
天 线 频谱 中 , 低频 的 峰值 大 约 比 高 频 的 平均 值 高 了 约 1 个 多 数量 级 , 而 在 两 天 线 互 相关 功率 
谱 中 ， 低 频 的 峰值 大 约 比 高 频 的 平均 值 高 了 约 2 个 多 数量 级 。 因 此 ， 从 互相 关 功 率 谱 出 发 ， 
更 有 利于 把 平台 干扰 信号 给 提取 出 来 。 


The Residual Data of A Antenna The Residual Data of B Antenna 
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Ca) A REG CLEAN 后 的 残余 信号 (b) B 天 线 CLEAN 后 的 残余 信号 
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图 11 A、B 天 线 CLEAN 后 的 残余 信号 


Fig. 11 The after-CLEAN residual signal of antennas A and B. (a)The residual signal of antenna A after CLEAN; (b) 
The residual signal of antenna B after CLEAN; 

我 们 进行 CLEAN 时 选取 的 6 增益 因子 为 0.2， 和 迭代 次 数 设置 为 10 万 次 。 从 经 过 CLEAN 之 
后 ， 天 线 A 和 天 线 了 的 残余 信号 显示 在 图 8 中 。 可 以 看 出 ，CLEAN 前 天 线 时 域 信号 中 存在 的 
一 些 强 峰 结构 已 经 消失 了 。CLEAN 之 后 天 线 A 数据 的 o 为 6.58 x 10-57， 天 线 B 数据 的 ao 为 
5.36 x 10V， 可 见 天 线 数 据 的 波动 性 已 经 得 到 极 大 地 减 小 。 天 线 信号 现在 在 不 同时 间 段 上 
的 幅度 近似 相等 ， 这 与 接收 机 的 白 噪声 本 底 的 特点 相符 合 。 


The model signal of A antenna The model signal of B antenna 
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0.4 | 
0.2 
0.2 
= = 
S 00 $ oo 
-0.2 
-0.2 
-0.4 
e | | | | -0.4 T T | 
T T T T T T T 
0 1000 2000 3000 4000 0 1000 2000 4000 
point(N) point(N) 
Ca) A 天 线 CLEAN 后 的 模型 信号 (b) B 天 线 CLEAN 后 的 模型 信号 


图 12 A、B 天 线 CLEAN 后 的 模型 信号 


Fig. 12 The model signal of antenna A and B after CLEAN. (a)Model signal of antenna A after CLEAN; (b) Model signal 


of antenna B after CLEAN; 
从 两 天 线 CLEAN Z Ja E SRA fa Ss CA 12) 中 可 以 看 出 ，CLEAN 之 后 的 model 
言 号 与 初始 信号 很 接近 ， 说 明 原始 的 信号 中 大 部 分 是 平台 干扰 信号 。 同 时 ， 两 天 线 model 
言 号 的 相关 系数 为 0.966， 接 近 于 1， 也 说 明 model 信号 中 大 部 分 是 相关 的 平台 干扰 信号 。 


The Spectrum of A Antenna After CLEAN The Spectrum of B Antenna After CLEAN 
| | 
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(a) A RBG CLEAN 后 频谱 Cb) B 天 线 CLEAN 后 频谱 


图 13 A、B 天 线 直接 CLEAN 后 频谱 


Fig. 13 Spectrum of antennas A and B after CLEAN. (a)Spectrum of antenna A after CLEAN; (b) Spectrum of antenna B 
after CLEAN; 

CLEAN 之 后 两 天 线 残 余 信 号 的 频谱 (图 13) 幅度 峰值 由 原来 的 约 1071V 降低 到 约 10-5V， 

降低 约 4 个 数量 级 。 残余 信号 的 频谱 也 不 再 有 明显 的 频谱 结构 , 在 各 个 频率 分 量 的 强度 大 致 

相等 ， 可 以 看 作 白 噪声 谱 。 这 说 明 ， 天 线 信号 中 的 强 平台 干扰 成 分 得 到 了 非常 有 效 的 去 除 。 
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The Cross-correlation spectrum After CLEAN 
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14 A. B KÉR CLEAN 后 的 来 态 互相 关 功 率 谱 
Fig. 14 The cross-power spectrum of antennas A and B after CLEAN 
CLEAN 后 的 互相 关 功 率 谱 《 图 14) 同样 印证 了 算法 的 有 效 性 。 可 以 看 出 ， 两 天 线 的 互相 
关 功 率 谱 中 已 经 没有 明显 的 峰值 结构 ， 各 个 频段 的 互相 关 功 率 谱 峰值 都 接近 于 10-?Y“， 对 
比 图 7 中 的 结果 ， 整 整 降低 了 约 9 个 数量 级 。 


the correlation coefficient of B and A antenna 


T T T T T T 
0 20000 40000 60000 80000 100000 
times(N) 


图 15 A、B 天 线 CLEAN 时 相关 系数 随和 迭代 次 数 分 布 曲线 

Fig.15 Distribution curve of correlation coefficient between antennas A and B vs. CLEAN iteration. 

我 们 给 出 了 两 天 线 信号 之 间 的 相关 系数 随 CLEAN 迭代 次 数 变 化 的 曲线 (图 15) o CLEAN 
之 前 初始 的 相关 系数 为 0.965， 由 强 平台 干扰 信号 所 主导 ; CLEAN 十 万 次 之 后 的 相关 系数 接 
近 于 0.2， 根 据 此 相关 系数 计算 出 的 信 噪 比 ， 可 以 反映 出 此 时 残余 信号 中 的 平台 干扰 信号 和 
接收 机 噪声 基本 在 同一 个 水 平 上 ， 此 时 信和 号 由 接收 机 噪声 所 主导 。 当 然 ，CLEAN 结束 之 后 ， 
两 天 线 信 号 的 相关 系数 并 没有 显著 接近 于 0。 这 说 明 ， 残 余 信号 中 可 能 包含 有 来 自 于 宇宙 的 
低频 射电 信号 ， 以 及 没有 彻底 清除 的 平台 干扰 信号 。 
4. 3 CLEAN 方法 效果 分 析 

我 们 假设 A、B 两 个 天 线 的 时 域 数据 分 别 表示 为 f4(t) = S(t) + N4(t) 和 f3(t) = S(t) + 
N3(t)， 其 中 S(t) 为 完全 相同 或 相关 的 信号 成 分 ， 设 其 均值 为 0%，，o 记 为 cs，N4(t) 和 N3(t) 分 
别 表示 A 和 B 天 线 的 接收 机 噪声 信号 ， 其 均值 也 都 为 0%，o 都 记 为 ov。 但 是 ，N4(t) 和 N3(t) 是 


互相 独立 的 。 相 关系 数 采用 线性 相关 系数 的 定义 r = ZED, [SW LESNRSE X gos/ou. BOSE 


推导 信 品 比 与 相关 系数 之 间 的 关系 ， 可 以 得 到 如 下 公式 : 
SNR? 
nium s 


当 相 关系 数 从 最 初 的 0. 965 降低 到 最 后 接近 0. 2 时 ， 由 公式 C 计算 得 到 的 信 噪 比 SNR 从 
5.3 降低 到 了 0.5。 说 明 在 末 态 残余 信号 中 ， 相 干 的 信号 成 分 和 不 相干 的 噪声 成 分 几乎 在 同 
信永 下 ES 
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通过 将 末 态 信号 转化 为 流量 密度 ， 得 到 B 天 线 末 态 信号 的 流量 密度 随 频率 的 分 布 〈 图 
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15) 。 从 中 可 见 ， 末 态 流量 密度 基本 上 降低 到 10-3W/(m? Hz 以下， 与 太阳 宁静 时 的 信和 号 


水 平 相当 。 后 续 可 以 在 此 基础 上 对 强 射电 爆发 


5 wit 


{Tt FE 


The Flux Density of antena B after CLEAN 


W/(m? *Hz2) 


r r T + 
1.0 10.75 20.5 30.25 
MHz 


3| 15 B 天 线 CLEAN 后 残余 信号 的 流量 密度 


Fig. 15 The residual flux density of antenna B after CLEAN 


经 过 CLEAN 算法 处 理 之 后 ,嫦娥 四 号 低频 射电 频谱 仪 各 天 线 的 探测 灵敏 度 ( 未 积分 ) 从 


10-5w/(m? - Hz) 提 高 到 10-3W/(m? Hz), BRS PMR. A 


E 此 灵敏 度 下 ， 可 以 对 太阳 


爆发 进行 监测 。 太 阳 爆 发 的 流量 密度 一 般 在 10-18W/(m2 - Hz)， 极 端 爆发 可 以 达到 10-1SW/ 


(m? - Hz)。 除 此 之 外 , 这 样 的 灵敏 度 也 足以 探测 部 分 强 的 快速 射 上 


He (Fast Radio Burst, FRB), 


其 辐射 流量 密度 可 以 达到 几 十 万 央 斯 基 (Jy) ，1HJy = 107-2°W/(m?- Hz。 如 果 再 对 时 域 数 


据 进行 积分 平均 ， 那 么 可 以 ; 


步 提 高 探测 灵敏 度 。 例 如 ， 如 果 对 一 道 4096 个 数据 i 


行 积 


分 平均 ， 则 灵敏 度 可 以 额外 提高 64 倍 。 这 将 大 大 拓展 低频 射电 频谱 仪 的 科学 探测 目标 。 
一 系列 用 CLEAN 算法 分 离 出 来 的 平台 干扰 信号 也 很 有 用 。 首 先是 被 用 来 分 析 平 台 干 扰 信 


号 的 频谱 、 


DAH 


空间 和 时 i 
磁 兼 容 设计 ， 大 幅度 降低 平台 干扰 水 平 。 
分 利用 这 些 知 识 ， 发 现 那 些 特征 异常 的 信号 ， 这 些 异常 信号 全 
上 干扰 信号 的 最 终 目的 是 进行 有 效 的 天 文 观测 。 那么 , 如 何在 强烈 的 平台 干扰 


其 次 ， 在 全 面 了 解 平 


司 变 化 特征 ， 并 确定 其 来 源 ， 以 便 在 今后 开发 类 似 设备 的 时 候 ， 改 进 电 


台 干 扰 信 号 特征 的 基础 上 , 充 


往来 自 于 宇宙 。 


背景 下 寻找 并 确 
扰 往往 有 固定 的 一 组 频谱 特征 ， 


的 射电 爆发 信号 。 如 果 这 些 频谱 还 具备 色散 等 特 和 


WS 


BASIE? ERIN Be Roa fu. qu. MERR. FAF 
旦 有 例外 的 频谱 信号 出 现 , 可 能 意味 着 记录 到 了 来 自 宇 宙 


E, 那么 几乎 可 以 肯定 它 来 


HFH. 再 例如 ， 


A、B、C 天 线 信 号 的 强度 比例 系数 。 若 平台 干扰 信号 产生 的 空间 分 布 是 固定 的 ， 那 么 不 管 其 


在 时 域 如 何 变化 ， 信 和 号 之 间 的 强度 比例 系数 也 是 


固定 的 。 倘 若 这 些 比例 系数 有 明显 的 改变 ， 


则 往往 意味 着 天 线 接收 到 了 来 自 宇宙 的 强 射 电 爆 发 信号 。 若 射电 爆发 信号 的 时 标 长 达 几 天 
(此 处 1 天 =24 小 时 ) 的 量 级 ， 那 么 这 些 比例 系数 还 会 随 着 月 球 自转 ， 发 生 规律 性 变化 ， 有 


助 于 更 加 可 靠 的 有 


的 数据 分 析 工 作 。 


AE x de 


6 结论 
本 论文 提 


的 方法 


相关 获得 可 见 度数 据 ， 对 可 见 度数 据 作 成 图 等 处 理 
进行 作 类 似 的 操作 ， 获 得 其 详细 的 频谱 、 


有 定 爆发 信号 的 来 源 。 当 然 ， 具体 如 何 搜 寻 异 党 


言 写 ， 还 得 依赖 于 后 续 细 致 


也 可 以 应 用 到 其 它 射 电 观测 设备 , 尤其 是 地 基 的 综合 射电 干涉 阵 。 H 


前 地 基干 涉 阵 受 到 越 来 越 多 的 来 自 地 面 设备 以 及 空间 卫星 的 无 线 电 干扰 。 借 助 CLEAN 算法 ， 
我 们 利用 俩 俩 天 线 之 间 的 互相 关 功 率 谱 , 在 信号 层面 把 干扰 信号 提取 出 来 , 再 对 残余 数据 做 
工作 。 同 样 ， 也 可 以 对 提取 出 的 干扰 信号 


空间 分 布 、 时 域 变化 等 信息 。 


一 种 方法 ， 它 从 两 组 信号 的 相关 性 出 发 ， 基 于 CLEAN 算法 ， 并 借助 互相 关 功 


率 谱 、 傅 立 叶 级 数 展开 等 工具 ， 把 低频 射电 频谱 仪 A、B、C 天 线 的 时 域 观测 数据 分 为 强 相关 


的 CLEAN 模型 信号 和 部 分 相关 的 残余 信号 两 个 部 分 。 其 中 CLEAN 模型 信号 主要 由 平台 干扰 信 
号 和 可 能 的 低频 强 射电 爆发 组 成 ; 残余 信号 由 接收 机 噪声 、 未 扣除 干净 的 平台 干扰 信号 和 常 
规 的 低频 射电 信号 组 成 。 将 该 方法 应 用 到 实际 数据 中 , 结果 表明 娣 娥 四 号 低频 射电 频谱 仪 的 
未 积分 灵敏 度 可 以 提高 约 8 个 数量 级 ， 达 到 10 7W/ mHz) 水平。 

在 此 基础 上 , 基于 对 平台 干扰 信号 中 确定 性 成 分 和 宽带 随机 成 分 的 分 别处 理 , 再 借助 低 
频 射 电 爆 发 信号 和 平台 干扰 信号 在 功率 谱 上 不 同 的 表现 , 以 及 常规 低频 射电 信号 受 月 球 自转 
调制 等 信息 , 将 来 工作 的 重点 是 进一步 分 析 CLEAN 模型 信号 和 残余 信号 ， 以 期 搜索 低频 强 射 
电 爆 发 信号 ， 乃 至 对 全 天 区 进行 粗略 的 成 像 。 
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Abstract: The Chang'e-4 low-frequency radio spectrometer is placed on the back of the 
moon, and its natural observation conditions are unique. However, the Chang'e-4 platform has 
strong interference of about 10 15*W/(m? - Hz) level, and there are obvious differences in the 
interference in each channel of time domain data, which greatly weakens the astronomical 
observation sensitivity of the low-frequency radio spectrum analyzer. To this end, starting from 
the correlation of the two sets of signals, this paper proposes to divide the time-domain 
observation data of the low-frequency radio spectrum analyzer A, B, and C antennas into strongly 
correlated CLEAN model signals and partial correlated residual signal based on the CLEAN 
algorithm, with the help of cross-correlation power spectrum, Fourier series expansion and other 
tools. The CLEAN model signal is mainly composed of platform interference signals and possible 
low-frequency strong radio bursts; the residual signal is composed of receiver noise, undeducted 
platform interference signals and conventional low-frequency radio signals. Applying this method 
to actual data, the results show that the unintegrated sensitivity of the Chang'e-4 low-frequency 
radio spectrum analyzer can be increased by about 8 orders of magnitude, reaching a level of 
10 75W/(m?- Hz). On this basis, based on the separate processing of the deterministic 
components and broadband random components in the platform interference signal, and then 
relying on the different performance of the low frequency radio burst signal and the platform 
interference signal in the power spectrum, and the conventional low frequency radio astronomy 
signal modulated by the moon rotation And other information, the focus of future scientific 
analysis work is to further process the CLEAN model signal and residual signal, with a view to 
discovering low-frequency strong radio astronomy burst signals, and even roughly imaging the 
entire sky. 
Key words: Chang'e-4 low-frequency radio spectrum analyzer; Platform interference removal; 
CLEAN algorithm 
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